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We investigate a mechanism of the appearance of kinetic friction in granular materials. We 
consider a small number of intervening inelastic particles between two rough surfaces as one of the 
simplest dynamical models to study granular friction. The resistance force applied to the upper 
surface is numerically calculated. We find that the resistance force F{t) is scaled as F'{vt) for a 
small pulling velocity v. The time average Fo — {F{t)} in the limit u ^ is not zero owing to 
the mutual collisions between the intervening particles. The nonzero Fq implies the appearance of 
kinetic friction in this simple dynamical system. 

PACS numbers: 



I. INTRODUCTION 

Friction forces have been long studied since the formulation of Coulomb- Amonton's laws. The laws of friction have 
been widely studied from the macroscopic scale to the nano scale J., 2, 3, 4]. The establishment of frictional contact 
between two solids by a normal force and its disruption by a tangential force are important processes in friction, and the 
frictional contact has been directly observed in experiments 'E?! . Recently, Miiser and coworkers and Daly et al. have 
proposed another mechanism of friction, that is, the origin of Coulomb- Amonton's laws is the presence of intervening 
dirt molecules or atoms between two solid surfaces @, 0, |8| . In a model of kinetic friction Q , Miiser assumed a potential 
between the intervening particle and the top (or bottom) surface such as Vt^t = Vt^b '^os{x/bt^b) + cos{2x /bt^b) + ■ ■ ■ , 
where 2-Kht.b are the periods of the top and bottom wall, respectively. A certain viscous force acting on the lubricant 
atom is further assumed for energy dissipation. A smooth motion is observed for the intervening particle and no 
kinetic friction appears, if the interaction force between the particle and the surface is a simple sinusoidal one, i.e., 
V("j = for n > 1; however, kinetic friction appears as a result of a stick-slip motion, if the interaction forces include 
a certain higher-harmonic component, e.g., V^f^ > 0, V/^ > 0, and Vj"^ = for n > 2. 

Recently, the friction and rheological properties of granular materials have been intensively studied. Granular 
materials exhibit various states such as a gaslike state, a liquidlike state and a solidlike state. They are numerically 
studied by molecular dynamics simulations including many inelastic particles. One of the interesting phenomena is the 
clustering in the gas phase owing to inelastic collisions [9]. Stick-slip motions caused by the frictional vibration have 
been studied numerically and experimentally p^ . [ill . . They are related to the transition between the liquidlike 
state and the solidlike state [l^, [3| . The transition from the liquidlike state to the solidlike state when the density is 
increased is called jamming transition p^ . Several scaling relations were found for the static properties such as the 
pressure and the average coordination number near the jamming transitio n Il6l. Scaling laws for rheological properties 
of sheared granular materials were also studied by numerical simulations [13, [3 • Below the critical density, the shear 
stress S is proportional to the square of shear velocity. This is called the Bagnold scaling [l^. Above the critical 
density, S is not zero even for an infinitesimal shear rate. The nonzero shear stress indicates the rigidity of the solid 
phase. A nonzero shear stress can be interpreted as the appearance of a kinetic friction in the granular materials. 
The nonzero shear stress appears even in an assembly of frictionless inelastic particles. That is, even if there is no 
friction among elemental particles, macroscopic friction appears. A force-chain network appears above the critical 
density in granular materials. It has been suggested that the anisotropy of the force-chain network is the origin of 
kinetic friction [20| . However, the origin of kinetic friction is not well understood, even in an assembly of such simple 
frictionless inelastic particles. A granular material begins to flow when the shear stress is beyond its critical value, 
called the yield stress in experiments, under a constant shear stress. The existence of yield stress is related to the 
appearance of the maximum static friction. 

In this study, we will mainly analyze kinetic friction by numerical simulation of a small number of inelastic inter- 
vening particles, expecting that some aspects of the friction may be clarified, if the degree of freedom of motion is 
reduced. Firstly, we neglect tangential forces such as the Mindlin force and rotational motion to simplify the system 
as much as possible. A small number of frictionless inelastic particles are intervening between the upper and lower 
walls. Our model is related to the Miiser model in that the intervening particles play an important role, but the 
model equation is that conventionally used in the numerical simulation of granular materials. The upper and lower 
walls are rough surfaces constructed from an array of semicircles. Inelastic collisions occur among the intervening 
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particles and between the intervening particles and the rough surfaces. We numerically calculate the resistance force 
exerted on the upper wall by the intervening particles, and discuss the origin of kinetic friction. In §2, we will discuss 
kinetic friction, static friction, and stick-slip motion in a single-particle system. In §3, we will study a two-particle 
system, and point out the importance of the mutual collisions between intervening particles for the appearance of 
kinetic friction. In actual experiments, the rotational degree of motion plays an important role [13, HH- In §4, we 
incorporate rotational motion by adding some tangential forces, and discuss several effects of rotational motion. In 
this study, no explicit form of friction is assumed between elemental particles. 

II. RESISTANCE FORCE BY A SINGLE INTERVENING PARTICLE 

We first consider a single-particle system confined between two walls. The particle is assumed to be a circle of radius 
r = ri. The upper and lower walls are constructed from periodic arrays of semicircles of radius r = r2 = 0.5 > ri as 
shown in Fig. 1(a). The y-coordinates of the centers of the semicircles in the lower and upper walls are respectively 
fixed to be and h > 0. We assume that the particle and semi-circles interact with a linear elastic force and a viscous 
force, only if the distance |xj — x^j between the particle i and the semicircle j is smaller than the sum of the radii of 
the two particles. Then, the intervening particle obeys the equation of motion; 

(fx- 

j 

where i — 1 indicates the particle intervening between the two walls, m is the mass of the particle, Sij = 1 if 
|xj — Xi| < rj + Ti, 5ij = if |xj — Xi| > Tj -t- ri, Wij = (xj — Xi)/|xj — x^j, and Vj — dxj/dt. A summation 
is taken for all the semicircles in contact with the particle. The reaction force exerted on the semicircles in the 
upper and lower walls by the particle is calculated to study friction force. The parameter values are ri = 0.425, 
TO = Ar\,K = 0.42 x 10^, and 77 = 83.9. A short-range viscous force effectively plays a role of an inelastic collision. 
Periodic boundary conditions are assumed at a; = 1 and a; = 2, that is, the position of the particle is reset to x — 1 
when X reaches the right boundary x = 2. Tangential forces can be assumed between two particles as elemental forces, 
which is more realistic for granular particles. However, we neglect such forces in this section for simplicity. 

Firstly, we performed a type of numerical simulation in which the upper wall is pulled at a constant velocity v in the 
x-direction, and the lower wall is fixed. The resistance force F exerted on the upper wall by the intervening particle 
was numerically evaluated. Figure 1(b) shows the time evolution of F atv = 0.02 for h = 0.56. The resistance force 
F{t) changes periodically with time. Figure 1(c) shows the time average {F) of the resistance force F{t) as a function 
of V for h — 1.68, 1.64, and 1.6. The average resistance force is roughly proportional to v for any h < he ^ 1.703. The 
intervening particle is trapped owing to the inelastic collisions by the lower wall, and the average velocity becomes 
for /i < he- This is partially related to the assumption that the rotational degree of freedom is neglected. The 
unidirectional motion induced by tangential forces is discussed later in §4. At the critical height, a semicircle in the 
upper wall is located just above the intervening particle, and the particle is located just between the two semicircles 
in the lower wall, as shown in Fig. 1(a). The centers of the two semicircles of radius r2 are denoted as A and C, and 
the center of the particle of radius ri is denoted as B; D is the contac t point of the t wo semicircles in the lower wall 
in Fig. 1(a). The critical height can be evaluated as he = ri + r2 + \/ {r\ -f- ri)^ — r| ^ 1.703, because AB= r\ -|- r2, 
BD= \/(ri -hri)2 -r^, and h^ =AD. 

We have also numerically calculated the normal force TV applied to the upper wall by the intervening particle. 
The time average of the normal force N increases as h decreases from he- Figure 1(d) shows ^ = {F)/{N) for 
h — 1.68, 1.64, and 1.6. The coefficients /i overlap rather well in a range of small v values for three different values of 
h, and fi approaches zero as w — *■ 0. Figure 1(e) shows the proportional constant = {F)/v for a small w as a function 
of h. The proportional constant ly increases as h decreases. 

The intervening particle is trapped between the two semicircles on the lower wall for h < he, and the upper wall 
slips on the particle. The average velocity of the intervening particle is therefore zero. From the Galilean symmetry, it 
is possible that the average velocity of the intervening particle is v and the lower wall slips below the moving particle. 
These two slip states are bistable for h < he- Figure 1(f) shows two trajectories of the particle during the slip process 
at = 0.005 and 3.5 for h = 1.56. At -y ~ 0, the trajectory is almost symmetric with respect to x = 1.5, and the 
average value of the resistance force F is almost zero from the symmetry. As v is increased, the x-coordinate of the 
trajectory shifts toward the right, then the average {F) becomes nonzero as a result of the asymmetric trajectory of 
the particle. Generally, the kinetic friction on a solid surface takes a nonzero value, even if v is sufficiently small. 
Thus, we conclude that no kinetic friction appears in this single-particle system. 

This single-particle system is a very simple dynamical system, but its dynamics is rather complex. If the interval h 
between the two walls is larger than the critical value h^, it is possible that the intervening particle is trapped in the 
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FIG. 1: (a) A particle is intervening between the upper and lower walls. The height h between the two walls is 1.703. (b) Time 
evolution of F{t) a.t v = 0.02 for h — 1.56. (c) Time average {F) of the resistance force as a function of v for h = 1.69, 1.64, 
and 1.59. (d) Ratio fi = {F)/{N) as a function of v for h — 1.69, 1.64, and 1.59. (e) Proportional constant u — {F)/v for a 
small w as a function of h. (f) Trajectories of the particle during the slip process at « = 0.01 (dashed curve) and 4 (solid curve) 
for h = 1.56. 
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FIG. 2: (a) Time evolution of z of the intervening particle at /i = 1.75. (b) Time evolution of F{t) a.t h — 1.75. (c) 
Double-logarithmic plot of time average {F) as a function of v. The dashed line denotes {F) oc v'^. 



lower wall (or the upper wall) and that the particle is never in contact with the upper wall (or the lower wall) . Then, 
the resistance force is completely zero. For general initial conditions, the intervening particle collides with the two 
walls in a chaotic manner. This might correspond to the phase of granular gas or granular liquid. Figure 2(a) shows 
the time evolution of the vertical coordinate z{t) of the intervening particle at h — 1.75. The intervening particle is 
located near z ~ 0.8 and z ^ 0.95 for most of the time. This corresponds to the two trapping states in the lower wall 
{z ~ 0.8) and upper wall (z ^ 0.95). The switching between the two trapping states occurs intermittently. Figure 
2(b) shows the time evolution of the resistance force F{t) at h = 1.75. The time evolution is also intermittent. The 
resistance force \F{t)\ takes a large value during the switching periods. Figure 2(c) shows a logarithmic plot of the 
long time average of F{t) as a function of v for h = 1.75. The solid line denotes (F) cx u^, which might correspond to 
the Bagnold scaling found in a numerical simulation of granular materials [Tsj . 

In most experiments on friction, a constant external force is applied to a solid body, and in some experiments, a 
solid body is pulled at a constant velocity via a spring. We can perform corresponding numerical simulations in our 
single-particle system. That is, the upper wall is pulled in the x-direction by a constant external force F^x- The 
upper wall moves as a rigid body. We have calculated the resistance force F{t) exerted on the upper wall by the 
intervening particle. If the external force is small, a stationary state is finally obtained. That is, the resistance force 
F is equal to Fg^ and the upper wall and the intervening particle do not move in the final stationary state. If Fg^ is 
larger than its critical value Fc, the upper wall begins to move toward the right. The critical force Fc is interpreted 
as the maximum static friction. The rhombi in Fig. 3(a) show the maximum static friction as a function of h. The 
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FIG. 3: (a) Maximum static friction Fc denoted by o in constant-force simulations, and the maximum value Fm denoted by 
+ of the time-periodic resistance force in constant-velocity simulations, (b) Coefficient of the maximum static friction as a 
function of h. 
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FIG. 4: (a) One-particle system with the upper wall coupled to a spring. The end point of the spring is pulled at a constant 
velocity v. (b) Time evolution of vt — x of the upper wall. A stick-slip motion appears. 



maximum static friction decreases monotonically with h. In numerical simulations of constant velocity, the resistance 
force F{t) changes periodically with time aX h < he as shown in Fig. 1(b). We calculated the maximum value Fm 
of the time-periodic resistance force F{t) as a function of h. The maximum value F„i for v — 0.02 is shown by -I- 
in Fig. 3(a). It is natural that the maximum static friction Fc is equal to Fm- Figure 3(b) shows the coefficient 
yU = Fc/N of the maximum static friction as a function of h. The coefficient /i is a slightly decreasing function of h. 
The coefficient of the maximum static friction is not zero even in this single-particle system. The maximum static 
friction corresponds to the existence of yield stress in granular materials. 

We also performed another type of numerical simulation. The upper wall is forced via a spring whose right end is 
pulled at a constant velocity w, as shown in Fig. 4(a). The center of mass of the upper wall obeys 

^~dt^ ^ ^ ^^^^ ~ ^^1) + ^^^^ ~ ' ■ + ^(^* - (2) 

id 

where X is the x-coordinate of the center of mass of the upper wall, is a unit vector in the x-direction, M is the 
mass of the upper wall, and k is the spring constant. Figure 4(b) shows the time evolution oi vt — X for M — 5, k — 5, 
and V = 0.08 at ft, = 1.5. A stick-slip motion is clearly seen. The upper wall almost sticks to the intervening particle 
and the lower wall for most of the time owing to the friction. The upper wall slips forward abruptly by AX ~ 21 in 
each period, which is much larger than the diameter of semicircles. The oscillation is not due to the periodic force 
given by the array of the semicircles, but is a friction-induced vibration observed in many frictional systems. Such 
frictional vibration occurs because the maximum static friction is larger than the kinetic friction, or the frictional force 
becomes weaker when the solid body begins to move. In our system, the kinetic friction is effectively zero; therefore, 
it is natural to observe the frictional vibration. This type of stick-slip motion is observed experimentally in granular 
systems. [13, Ell 



III. KINETIC FRICTION BY INTERVENING PARTICLES 



As shown in the previous section, no kinetic friction appears in a single-particle system. We consider a two-particle 
system in this section. To investigate kinetic friction, the upper wall is pulled at a constant velocity v. Two particles 
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FIG. 5: (a) Resistance force {F) as a function of v for h — 1.56, 1.51, and 1.46. (b) Kinetic friction force Fo for infinitesimally 
small V as a function of h. (c) The coefficient of the kinetic friction as a function of h 




FIG. 6: (a) Time evolutions of and X2{t) for h — 1.46 and v — 0.02. (b) Time evolution of d = 

^y{xl{t) - X2{t)y + {yi{t) - y2{t)y ~ 2ri for h = 1.46 and v = 0.02. (c) Snapshot patterns for h = 1.46 and v = 0.02 at 
t = 2400 (top) and t = 2410 (bottom). 



with the same radius r = ri — 0.425 are intervening between the upper and lower walls. The periodic boundary 
conditions are imposed at a; = 1 and x = 3 for this two-particle system. In the two-particle system, the intervening 
particles interact with each other with the same types of elastic and viscous forces, only when the distance between 
the two particles is smaller than 2ri. The two intervening particles also interact with the upper and lower walls. The 
initial conditions are xi — 1.05, j/i = h/2,Vxi — Vyi — —0.05,2:2 — 2.102, ?;2 ~ h/2, and Vx2 ~ Vy2 = —0.05 for the 
two intervening particles. The initial conditions for the semicircles in the lower wall are Xj = j^yj = 0, and Vxj = 0. 
For the semicircles in the upper wall, Xj = j + 0.2505, yj = h, and Vxj = v = 0.01. The velocity v of the upper wall is 
increased stepwise to 0.01 x (n + 1) at t = 2500 x n (n = 1, 2, • • • ). There is another critical width hc2 ^ 1-55 in this 
two-particle system. For hc2 < h < he, the average resistance force (F) per particle is roughly proportional to u, which 
is similar to the behavior of the single-particle system. This is because the two particles are almost trapped between 
different neighboring semicircles in the lower wall, and there is no interaction between the two intervening particles. 
However, the behavior of the resistance force changes when h is smaller than hc2- Figure 5(a) shows the average 
resistance force (F) as a function of v for h — 1.56, 1.51, and 1.46 for the two-particle system. At h — 1.56 > hc2, the 
average resistance force is roughly proportional to the velocity v. At h — 1.51 and 1.46, the average resistance force 
takes a nonzero value, Fq, even for a sufficiently small v. The nonzero value Fq in the limit w — > appears when h 
is smaller than hc2 ^ 1.55. We can interpret this resistance force as a kind of kinetic friction. Figure 5(b) shows _fo 
as a function of h. We have calculated the normal force acting on the upper wall and evaluated the kinetic friction 
coefficient fig for a sufficiently small v. Figure 5(c) shows the coefficient /i^ as a function of h. On the other hand, the 
coefficient of the maximum static friction for the two-particle system was the same as that in Fig. 3(b). The kinetic 
friction coefficient is definitely smaller than the coefficient of the maximum static friction in our system. 

The appearance of a nonzero value Fq is closely related to the mutual collisions of the two intervening parti- 
cles. Figure 6(a) shows the time evolutions of Xi{t) and X2{t) for h = 1.46 and v = 0.02, where Xi (solid curve) 
and X2 (dashed curve) denote the x-coordinates for the first and second particles correspondingly; however, no re- 
setting process x — > x — 2 under the periodic boundary conditions is taken into account in the notation x. The 
dotted line represents x ~ 0.01333i. The intervening particles are moving at an average velocity 0.0133, which is 
different from the velocities of 0.02 and of the upper and lower walls. Figure 6(b) shows the time evolution of 
d = ^y {xi{t) — X2(t))'^ + {yi(t) — y2{t)Y ~ 2ri. Mutual collision occurs when d is negative. Figure 6(c) shows two 
snapshots of the two particles at t = 2400 and 2410. The mutual collision is observed in the snapshot at i = 2410. 

Figure 7(a) shows the average velocity V of xi as a function of v. The average velocity satisfies V = (2/3)w. The 
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FIG. 7: (a) Average velocity V of ii as a function of v. (b) Time evolution of the resistance force F at v — 0.03 (solid curve) 
and F' — F{t/3 + 200/3) at u = 0.01 (dashed curve). The difference is hardly observed. 




FIG. 8: (a) Snapshot pattern of a ten-particle system at t = 2500 for v = 0.03 and h — 1.3. (b) Time evolution of Xi for 
i = 1, 2, ■ ■ ■ , 10 at u = 0.03 and h — 1.3. The dashed line denotes x ~ {3/7)vt. (c) Spatio-temporal plot of mutual collisions. 
Neighboring particles satisfying d — \J (x\{t') — x-iit))^ + {yi{t) — y2(t)^ ~ 2ri < at time t are linked with a line. 



ratio of the velocity of the intervening particles to that of the upper wall is a simple fraction, which implies a type of 
locking phenomenon. The locking phenomenon appears because of the spatially periodic structure of the upper and 
lower walls. These locking behaviors depend on the initial conditions. That is, there are multistable states in this 
system. The trapping states with V — or V — v are also stable states even for h < hc2- We have confirmed the 
existence of other locking states with V — {l/2)v and V — (5/6)w for the same h. 

To understand the appearance of the nonzero value Fq for a sufficiently small v, we have investigated the time 
evolution of the resistance force F{t) applied to the upper wall in more detail. The solid curve in Fig. 7(b) shows 
the time evolution of F{t) at v = 0.03. The resistance force changes periodically with time, and takes positive and 
negative values. The behavior is similar to that in the case of the single-particle system. However, at the moment of 
mutual collisions, impulsive forces are exerted on the lower and upper walls by the intervening particles. In the time 
evolution of F{t) in Fig. 7(b), the impulsive forces appear as the second peak in the doublet peaks in the range of 
F > 0. As a result of the additional impulses, the time average of F{t) takes a positive value (F) = 11.7 at v = 0.03. 
The dashed curve in Fig. 7(b) shows the time evolution of the resistance force at t; = 0.01. However, the time scale 
is reduced by one-third. That is, F' = F{t/3 + 200/3) is shown in this figure. The solid and dashed curves are well 
overlapped. This implies that the time evolution is scaled as F{t) = F'{vt) for a sufficiently small v, and therefore, 
the time average (F) of F is nonzero even if v is sufficiently small. That is, a certain number of impulsive forces are 
applied to the upper wall, while the upper wall shifts by one period in the x-direction. For a sufficiently small v, the 
process occurs slowly, but the total force applied to the upper wall per spatial period does not depend on the velocity 
V. This is the origin of the appearance of a nonzero Fq in our two-particle system. 

Furthermore, we have studied kinetic friction in some many-particle systems. Figure 8 shows some numerical results 
for a ten-particle system at h — 1.3. Figure 8(a) displays a snapshot pattern of the ten particles at t — 2500 for 
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FIG. 9: (a) Average resistance force (F) for a single-particle system as a function of the velocity i; of the upper wall at 
h = 1.68, 1.64, and 1.6 for rjt = O.I77, and Kt — 0. (b) Average velocity F as a function of the velocity v of the upper wall at 
h = 1.68, 1.64, and 1.6. (c) Average velocity V of the intervening particle as a function of e = rjt/rj at v — 0.05 for h — 1.68 
and 1.6. The parameter e is increased stepwise by Ae = 0.0005 from e = 0. 



V = 0.03. Three intervening particles located near x 2 are in contact with each other in the snapshot. Figure 
8(b) shows the time evolution of for ? = 1, 2, • ■ • , 10. All intervening particles move at the same average velocity 

V ~ (3/7)f . This is another example of the locking phenomenon. Figure 8(c) shows a spatio-temporal plot of mutual 
collisions. That is, two particles satisfying d = y/ {xi{t) — X2(i))^ + {yi{t) — 2/2(0^ ~ < are linked by a line at 
time t. Clusters of two, three, or four intervening particles are created and annihilated with time by mutual collisions. 
The mutual collisions occur sequentially, and they propagate in the x-direction. The nonzero resistance force Fq is 
also observed in this system. The impulsive forces induced by the mutual collisions between intervening particles 
might be the origin of kinetic friction. 

IV. KINETIC FRICTION BY INTERVENING PARTICLES WITH ROTATIONAL DEGREES OF 

FREEDOM 

In the previous sections, we have neglected tangential forces. If the tangential forces are neglected, no torque is 
applied to each particle and then the rotational degree of freedom can be neglected. In this section, we incorporate a 
tangential elastic force and a tangential viscous force. The equation of motion of the intervening particle is written as 

(fx 

j 

+{KtA{{rj - Ti) ■ Uj + Tidi + r.j9j) + ?;t((vj - v^) • -I- r^w^ + rjUjj)}tij\5ij, 

(3) 

where Kt is the elastic constant for the tangential displacement, Oi is the rotation angle of the ith particle, A((rj — r^) ■ 
tij +ri9i +rjdj) denotes the tangential displacement from the instance of contact between two particles, oji = d6i/dt, 
5ij = 1 if |xj — Xi| < rj -\- ri, Sij = if |xj — x^j > rj + r^, and is a unit vector perpendicular to n^^. The angle 
9i and uji for the semicircles in the upper and lower walls are set to be 0, because the lower wall is stationary and 
the upper wall is shifted at a constant velocity. If the rotational degree of freedom is taken into consideration, the 
equation of rotation is written as 

^ ^ {^t^((rj - r,) • ty + + rj6j) + 77t((vj - v,) ■ ty + nuji + rjUjj)}5ij, (4) 

where / is the inertial moment of the ith particle. The inertial moment is assumed to be / = {l/2)mrf in this 
simulation. Firstly, we neglect the tangential elastic force, that is, Kt = 0. Figure 9(a) shows the relation of the 
average resistance force (F) for a single-particle system and the velocity v of the upper wall at /i = 1.68, 1.64, and 
1.6 for T]t — O.I77. Figure 9(a) corresponds to Fig. 1(c) for rjt — 0; however, the average resistance force F is larger 
than that in the case of rjt — 0, probably because of the additional viscous force. However, the resistance force (F) 
becomes zero when v approaches 0. That is, no kinetic friction appears in the single-particle system, even if rotational 
motion is taken into consideration. Figure 9(b) shows the average velocity V of the intervening particle as a function 
of the velocity v of the upper wall. The average velocity V satisfies V = v/2 for j]t — O.lry. For the case of rjt — 0, the 
velocity of the intervening particle is zero, that is, the intervening particle is trapped in the lower wall. The tangential 
viscous force induces the flow of the intervening particle. Figure 9(c) shows the average velocity V of the intervening 
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FIG. 10: (a) Average resistance force (F) for a two-particle system as a function of the velocity v of the upper wall at 
h = 1.56, 1.51, and 1.46 for rjt ~ O.lry. (b) Average velocity V of the first particle as a function of the velocity v of the upper 
wall a.t h — 1.56, 1.51, and 1.46 for rji — O.lr]. 
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FIG. 11: (a) Average resistance force (F) for a single-particle system as a function of the velocity i; at ft = 1.64 and 1.68 
for Kt — O.OIK and rjt — rj. The rotational degree of freedom is not taken into consideration, (b) The same as (a), but the 
rotational degree of freedom is taken into consideration. 



particle as a function oi e — rjt /r] a.t v — 0.05 for h — 1.68 and 1.6, when e is increased stepwise by Ae — 0.0005 
from 0. A discontinuous transition from the trapped state to the flowing state is observed. The critical values Cc are 
0.056 for h = 1.68, 0.0535 for h = 1.64, and 0.05 for h = 1.6. When e is decreased stepwise by Ae = 0.0005 from 0.1, 
another discontinuous transition is observed. The critical values ec2 are 0.0555 for h = 1.68, 0.053 for h = 1.64, and 
0.035 for h = 1.6. The transitions are almost discontinuous for h = 1.68 and 1.64, but a hysteresis is observed for 
h — 1.6. Figure 10 shows numerical results for a two-particle system at rjt — O.lr/. Figure 10(a) shows the average 
resistance force as a function of v for h = 1.46, 1.51, and 1.56. Figure 10(b) shows the average velocity V of the 
first intervening particle as a function of w for ft, = 1.46,1.51, and 1.56. The average resistance force (F) becomes 
zero for Tj ^ at ft = 1.56. The average velocity is = v/2 for v < 0.34. At ft = 1.56, no mutual collisions between 
the intervening particles occur, and each particle behaves independently. However, a nonzero average resistance force 
appears for ft — 1.46 and 1.51. This implies the existence of the kinetic friction induced by mutual collisions between 
the intervening particles. The average velocity also satisfies V — 2v/3 for a small v. The relation V — 2v/3 is not 
satisfied for a relatively large u at ft = 1.46 and ft = 1.51, but other ratios such as V/v — 1, 5/6, and 1/3 are observed. 
It might be an origin of the non-monotonic relation of (F) and ft. 

Finally, we consider the effect of the tangential elastic force for a single-particle system. First, the rotational degree 
of freedom is not taken into consideration in the numerical simulation, that is, 9i is set to be zero at any t. Figure 
11(a) shows the average resistance force (F) as a function of z; at ft = 1.68 and 1.64 for Kt = O.OIK and rjt = rj. The 
appearance of a nozero Fq in the limit u — *■ is seen in this figure. The nonzero Fq increases with Kt- The tangential 
elastic force seems to be an origin of kinetic friction in this single-particle system. However, if the rotational degree of 
freedom is taken into consideration, that is, 9i changes with time according to Eq. (4), the average resistance force Fq 
becomes zero for a sufficiently small v. Figure 11(b) shows the numerical results at ft = 1.68 and 1.64 for Kt = O.OIK 
and Tjt = rj. Thus, we conclude again in our numerical simulation that no kinetic friction appears in a single-particle 
system. 
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V. SUMMARY AND DISCUSSION 



We have studied frictional forces using a few inelastic particles intervening between two rough surfaces. The simplest 
model is a singlc-particlc system without the rotational degree of freedom. Even in this simplest system, the dynamical 
behavior is rather complex. We have found a kind of jamming transition by changing the interval h between the two 
walls. Below the critical interval, the intervening particle is trapped in the lower (or upper) wall. Time-averaged 
resistance force {F) is roughly proportional to the constant velocity v of the upper wall. No kinetic friction appears in 
a single-particle system, because {F) approaches zero in the limit v — » 0. At /i > /ic, a chaotic time evolution is found, 
and the average resistance force is approximated as {F) oc v^. We have confirmed that the maximum static friction 
is equal to the maximum value of the time-periodic resistance force F{t). The maximum static friction corresponds 
to the yield stress for the granular solidlike state. We have also confirmed the appearance of a stick-slip motion even 
in this single-particle system, when the upper wall is forced via a spring. 

Kinetic friction has appeared in a two-particle system. The mutual collisions between the intervening particles 
seem to play an important role in the appearance of kinetic friction and the nontrivial average velocity V = {2/3)v. 
At the moment of mutual collisions, impulsive forces are exerted on the lower and upper walls. The amplitude of the 
impulsive forces does not decrease even for a sufficiently small v. The resistance force F{t) behaves as F{t) = F'{vt) 
for a sufficiently small v; therefore, the time average of F{t) takes a nonzero constant value even for a sufficiently 
small V. We have further observed the propagation of mutual collisions in a ten-particle system. The mutual collisions 
between the intervening particles occur even for a sufficiently small pulling velocity. The mutual collisions between 
intervening particles induce the dissipation of kinetic energy in the x direction, which might lead to kinetic friction; 
however, the detailed mechanism of this phenomenon is not yet completely understood theoretically. 

Our model is similar to Miiser's model in that the intervening particles play an important role. However, no kinetic 
friction appears if there is only one intervening particle, but kinetic friction appears as a result of the mutual inelastic 
collisions between the intervening particles. We consider that the origin of kinetic friction in our model is different 
from that in Miiser's model. 

Finally, we have incorporated tangential viscous and elastic forces in the model equations. If the tangential viscous 
forces are included, the intervening particle exhibits a unidirectional motion with V = (l/2)v, when the tangential 
viscosity is larger than its critical value. However, the appearance or nonappearance of kinetic friction do not depend 
on the tangential viscous force. If both the tangential viscous force and tangential elastic force are included and the 
rotational degree of freedom is not taken into consideration, the nonzero resistance force Fq appears. However, if the 
rotational degree of freedom is taken into consideration, (F) becomes zero for a sufficiently small v. 

Our model system is simple, but its dynamical behavior is rather complex. It is left for future study to analyze the 
system from the viewpoint of the theory of dynamical system, as well as to associate the present results for a small 
number of particles with the dynamical behavior of many-particle systems. 
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